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Abstract 

We present a measurement of the W boson mass using data collected by 
the D0 experiment at the Fermilab Tevatron during 1994-1995. We identify 
W bosons by their decays to eu final states where the electron is detected 
in a forward calorimeter. We extract the W boson mass, M^r, by fitting 
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the transverse mass and transverse electron and neutrino momentum spectra 
from a sample of 11,089 W ^ ev decay candidates. We use a sample of 
1,687 dielectron events, mostly due to Z — > ee decays, to constrain our model 
of the detector response. Using the forward calorimeter data, we measure 
Mw = 80.691 ±0.227 GeV. Combining the forward calorimeter measurements 
with our previously published central calorimeter results, we obtain M\y = 
80.482 lb 0.091 GeV. 
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I. INTRODUCTION 

In this article we describe the first measurement of the mass of the W boson using 
electrons detected at large rapidities (i.e. between 1.5 and 2.5). We use data collected in 
1994-1995 with the D0 detector [§] at the Fermilab Tevatron pp collider. This measurement 
performed with the D0 forward calorimeters complements our previous measurements 
with central electrons and the more complete combined rapidity coverage gives useful 
constraints on model parameters that permit reduction of the systematic error, in addition 
to increasing the statistical precision. 

The study of the properties of the W boson began in 1983 with its discovery by the UAl 
and UA2 collaborations at the CERN pp collider. Together with the discovery of the 
Z boson in the same year [§,0, it provided a direct confirmation of the unified description 
of the weak and electromagnetic interactions |]10[, which — together with the theory of 
the strong interaction, quantum chromodynamics (QCD) — now constitutes the standard 
model. 

Since the W and Z bosons are carriers of the weak force, their properties are intimately 
coupled to the structure of the model. The properties of the Z boson have been studied in 



great detail in e^e~ collisions ||TT|. The study of the W boson has proven to be significantly 
more difficult, since it is charged and so cannot be resonantly produced in e~^e~ collisions. 
Until recently its direct study has therefore been the realm of experiments at pp colliders 
PJ5| JT^JTB[ ] . Direct measurements of the W boson mass have also been carried out at LEP2 
||1^^17| using nonresonant W pair production. A summary of these measurements can be 



found in Table XI at the end of this article. 



The standard model links the W boson mass to other parameters, 

^ / vra(M|) \ Ml 

[ V2Gf J (M| - M^){1 - /^TEw) ^ ' 

in the "on shell" scheme |T^. Aside from the radiative corrections ^tewi the W boson 
mass is thus determined by three precisely measured quantities, the mass of the Z boson 



Mz [0], the Fermi constant Gp [[T^j ^-^d the electromagnetic coupling constant a evaluated 
at Q2 = Ml lH: 

Mz = 91.1867 ± 0.0021 GeV, (2) 
Of = (1.16639 ± 0.00001) x 10"^ GeV'^ , (3) 
a = (128.88 ±0.09)-^ . (4) 

From the measured W boson mass, we can derive the size of the radiative corrections Arpw- 
Within the framework of the standard model, these corrections are dominated by loops 
involving the top quark and the Higgs boson (see Fig. 0). The correction from the tb loop is 
substantial because of the large mass difference between the two quarks. It is proportional 
to mf for large values of the top quark mass rrit. Since rrit has been measured pO| , pT| , this 
contribution can be calculated within the standard model. For a large Higgs boson mass, 
tuh, the correction from the Higgs loop is proportional to ln{mH)- In extensions to the 
standard model, new particles may give rise to additional corrections to the value of Mw- 
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FIG. 1. Loop diagrams contributing to the W boson mass. 

In the minimal supersymmetric extension of the standard model (MSSM), for example, 
additional corrections can increase the predicted W mass by up to 250 MeV p2| . 

A measurement of the W boson mass therefore constitutes a test of the standard model. 
In conjunction with a measurement of the top quark mass, the standard model predicts M^y 
up to a 200 MeV uncertainty due to the unknown Higgs boson mass. By comparing the 
standard model calculation to the measured value of the W boson mass, we can constrain the 
mass of the Higgs boson, the agent of the electroweak symmetry breaking in the standard 
model that has up to now eluded experimental detection. A discrepancy with the range 
allowed by the standard model could indicate new physics. The experimental challenge is 
thus to measure the W boson mass to sufficient precision, about 0.1%, to be sensitive to 
these corrections. 

II. OVERVIEW 
A. Conventions 

We use a Cartesian coordinate system with the z-axis defined by the direction of the 
proton beam, the x-axis pointing radially out of the Tevatron ring, and the y-axis pointing 
up. A vector p is then defined in terms of its projections on these three axes, p^^, Py, p^. 
Since protons and antiprotons in the Tevatron are unpolarized, all physical processes are 
invariant with respect to rotations around the beam direction. It is therefore convenient to 
use a cylindrical coordinate system, in which the same vector is given by the magnitude of 
its component transverse to the beam direction, p^, its azimuth (p, and p^- In pp collisions, 
the center-of-mass frame of the parton-parton collisions is approximately at rest in the plane 
transverse to the beam direction but has an undetermined motion along the beam direction. 
Therefore the plane transverse to the beam direction is of special importance, and sometimes 
we work with two-dimensional vectors defined in the x-y plane. They are written with a 
subscript T, e.g. pt- We also use spherical coordinates by replacing pz with the polar angle 
9 (as measured between p^ and the z-axis) or the pseudorapidity r] = — Intan (6'/2). The 
origin of the coordinate system is in general the reconstructed position of the pp interaction 
when describing the interaction, and the geometrical center of the detector when describing 
the detector. For convenience, we use units in which c = h = 1. 

B. W and Z Boson Production and Decay 

In pp collisions at ^/s = 1.8 TeV, W and Z bosons are produced predominantly through 
quark- ant iquark annihilation. Figure ^ shows the lowest-order diagrams. The quarks in the 
initial state may radiate gluons which are usually very soft but may sometimes be energetic 
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enough to give rise to hadron jets in the detector. In the reaction, the initial proton and 
antiproton break up and the fragments hadronize. We refer to everything except the vector 
boson and its decay products collectively as the underlying event. Since the initial proton 
and antiproton momentum vectors add to zero, the same must be true for the vector sum 
of all final state momenta and therefore the vector boson recoils against all particles in the 
underlying event. The sum of the transverse momenta of the recoiling particles must balance 
the transverse momentum of the boson, which is typically small compared to its mass but 
has a long tail to large values. 




u/d 




d/ u/d 

FIG. 2. Lowest order diagrams for W and Z boson production. 



We identify W and Z bosons by their leptonic decays. The D0 detector (Sec. p|) is 
best suited for a precision measurement of electrons and positronsQ and we therefore use 
the decay channel W ^ eu to measure the W boson mass. Z ^ ee decays serve as an 
important calibration sample. About 11% of the W bosons decay to eu and about 3.3% of 
the Z bosons decay to ee. The leptons typically have transverse momenta of about half the 
mass of the decaying boson and are well isolated from other large energy deposits in the 
calorimeter. Gauge vector boson decays are the dominant source of isolated high-p^ leptons 
at the Tevatron, and therefore these decays allow us to select clean samples of W and Z 
boson decays. 



C. Event Characteristics 

In events due to the process pp —>■ {W —>■ ev) + X, where X stands for the underlying 
event, we detect the electron and all particles recoiling against the W boson with pseudora- 
pidity —4 < r] < A. The neutrino escapes undetected. In the calorimeter we cannot resolve 
individual recoil particles, but we measure their energies summed over detector segments. 
Recoil particles with \rj\ > 4 escape unmeasured through the beampipe, possibly carrying 
away substantial momentum along the beam direction. This means that we cannot measure 
the sum of the z-components of the recoil momenta, Uz-, precisely. Since these particles es- 
cape at a very small angle with respect to the beam, their transverse momenta are typically 
small and neglecting them in the sum of the transverse recoil momenta, ut causes a small 
amount of smearing of ut- We measure ut by summing the observed energy flow vectorially 



In the following we use "electron" generically for both electrons and positrons. 
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over all detector segments. Thus, we reduce the reconstruction of every candidate event to 
a measurement of the electron momentum p{e) and ut- 

Since the neutrino escapes undetected, the sum of all measured final state transverse 
momenta does not add to zero. The missing transverse momentum ^t, required to balance 
the transverse momentum sum, is a measure of the transverse momentum of the neutrino. 
The neutrino momentum component along the beam direction cannot be determined, be- 
cause Uz is not measured well. The signature of a. W ^ eu decay is therefore an isolated 
high-pr electron and large missing transverse momentum. 

In the case of Z ^ ee decays, the signature consists of two isolated high-p-r electrons 
and we measure the momenta of both leptons, p(ei) and p(e2), and ut in the detector. 



D. Mass Measurement Strategy 

Since Pz{t^) is unknown, we cannot reconstruct the ez/ invariant mass for W ^ ev candi- 
date events and therefore must resort to other kinematic variables for the mass measurement. 
For recent measurements |12jT3|j5|,^ the transverse mass. 



rriT 



^2pT(eKM(l -cos (0(e) -0(z/))) , (5) 



was used. This variable has the advantage that its spectrum is relatively insensitive to the 
production dynamics of the W boson. Corrections to ttlt due to the motion of the W are of 
order (g^/M^y)^, where qx is the transverse momentum of the W boson. It is also insensitive 
to selection biases that prefer certain event topologies (Sec. |V1 Dp . However, it makes use 
of the inferred neutrino pt and is therefore sensitive to the response of the detector to the 
recoil particles. 

The electron pt spectrum provides an alternative measurement of the W mass. It is 
measured with better resolution than the neutrino pt and is insensitive to the recoil mo- 
mentum measurement. However, its shape is sensitive to the motion of the W boson and 
receives corrections of order qt/Mw- It thus requires a better understanding of the W boson 
production dynamics than the rriT spectrum does. 

These effects are illustrated in Figs. ^ and which show the effect of the motion of 
the W bosons and the detector resolutions on the shapes of the m-p and pr^e) spectra. 
The solid line shows the shape of the distribution before the detector simulation and with 
gT=0. The points show the shape after qx is added to the system, and the shaded histogram 
also includes the detector simulation. We observe that the shape of the rriT spectrum is 
dominated by detector resolutions and the shape of the Pt{g) spectrum by the motion of 
the W boson. 

The shape of the neutrino pt spectrum is sensitive to both the W boson production 
dynamics and the recoil momentum measurement. By performing the measurement using 
all three spectra, we provide a powerful cross check with complementary systematics. 

All three spectra are equally sensitive to the electron energy response of the detector. We 
calibrate this response by forcing the observed dielectron mass peak in the Z ^ ee sample 
to agree with the known Z mass |Tl| (Sec. |VT|). This means that we effectively measure the 
ratio of W and Z masses, which is equivalent to a measurement of the W mass because the 
Z mass is known precisely. 
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FIG. 3. The uit spectrum for W bosons with qt 
(•), and with detector resolutions (shaded). 
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To carry out these measurements, we perform a maximum likelihood fit to the spectra. 
Since the shape of the spectra, including all the experimental effects, cannot be computed 
analytically, we need a Monte Carlo simulation program that can predict the shape of the 
spectra as a function of the W mass. To measure the W mass to a precision of order 100 MeV, 
we wish to estimate individual systematic effects with a statistical error of 5 MeV. Our 
technique requires a Monte Carlo sample of 10 million accepted W bosons for each such effect. 
The program therefore must be capable of generating large event samples in a reasonable 
time. We obtain the required Monte Carlo statistics by employing a parameterized model 
of the detector response. 

We next summarize the aspects of the accelerator and detector that are important for our 
measurement (Sec. |ITTp. Then we describe the data selection (Sec. |Vp and the fast Monte 
Carlo model (Sec. 0). Most parameters in the model are determined from our data. We 
describe the determination of the various components of the Monte Carlo model in Sees. [VJ- 
IX| . After tuning the model, we fit the kinematic spectra (Sec. ^), perform some consistency 
checks (Sec. ^), and discuss the systematic uncertainties (Sec. pCIip . We present the error 
analysis in Sec. pCIII|, and summarize the results and present the conclusions in Sec. pCIV. 



III. EXPERIMENTAL METHOD 
A. Accelerator 



During the data run, the Fermilab Tevatron collided proton and antiproton beams 
at a center-of-mass energy of ^/s = 1.8 TeV. Six bunches each of protons and antiprotons 
circulated around the ring in opposite directions. Bunches crossed at the intersection re- 
gions every 3.5 fis. During the 1994-1995 running period, the accelerator reached a peak 
luminosity of 2.5 x 10^^cm~^s^^ and delivered an integrated luminosity of about 100 pb~^. 
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FIG. 4. The pT{e) spectrum for W bosons with qx = ( ), with the correct qj- distribution 

(•), and with detector resolutions (shaded). 

The beam interaction region at D0 was at the center of the detector with an r.m.s. length 
of 27 cm. 

The Tevatron tunnel also housed a 150 GeV proton synchrotron, called the Main Ring, 
used as an injector for the Tevatron and accelerated protons for antiproton production 
during collider operation. Since the Main Ring beampipe passed through the outer section 
of the D0 calorimeter, passing proton bunches gave rise to backgrounds in the detector. We 
eliminated this background using timing cuts based on the accelerator clock signal. 

B. Detector 

1. Overview 

The D0 detector consists of three major subsystems: an inner tracking detector, a 
calorimeter, and a muon spectrometer. It is described in detail in Ref. [0]. We describe only 
the features that are most important for this measurement. 

2. Inner Tracking Detector 

The inner tracking detector is designed to measure the trajectories of charged particles. 
It consists of a vertex drift chamber, a transition radiation detector, a central drift chamber 
(CDC), and two forward drift chambers (FDC). There is no central magnetic field. The 
CDC covers the region \ri\ < 1.0. The FDC covers the region 1.4 < |?7| < 3.0. Each FDC 
consists of three separate chambers: a $ module, with radial wires which measures the 
coordinate, sandwiched between a pair of G modules which measure (approximately) the 
radial coordinate. Figure ^ shows one of the two FDC detectors. 
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3. Calorimeter 



The uranium/liquid-argon sampling calorimeter (Fig. ^ is the most important part of 
the detector for this measurement. There are three calorimeters: a central calorimeter (CC) 
and two end calorimeters (EC), each housed in its own cryostat. Each is segmented into 
an electromagnetic (EM) section, a fine hadronic (FH) section, and a coarse hadronic (CH) 
section, with increasingly coarser sampling. 




FIG. 6. A cutaway view of the D0 calorimeter and tracking system. 



12 



The ECEM section (Fig. 0) has a monohthic construction of alternating uranium plates, 
liquid-argon gaps, and multilayer printed-circuit readout boards. Each end calorimeter is 
divided into about 1000 pseudo-projective towers, each covering 0.1x0.1 in r/ x 0. The EM 
section is segmented into four layers, 0.3, 2.6, 7.9, and 9.3 radiation lengths thick. The 
third layer, in which electromagnetic showers typically reach their maximum, is transversely 
segmented into cells covering 0.05x0.05 in rj x (p. The EC hadronic section is segmented 
into five layers. The entire calorimeter is 7-9 nuclear interaction lengths thick. There are no 
projective cracks in the calorimeter and it provides hermetic and almost uniform coverage 
for particles with |?7| < 4. 




D0 END CALORIMETER ELECTROMAGNETIC MODULE 

FIG. 7. The ECEM section of an end calorimeter. 

The signals from arrays of 2x2 calorimeter towers covering 0.2x0.2 in rj x (p are added 
together electronically for the EM section alone and for the EM and hadronic sections 
together, and shaped with a fast rise time for use in the Level 1 trigger. We refer to these 
arrays of 2x2 calorimeter towers as "trigger towers." 

The liquid argon has unit gain and the end calorimeter response was extremely stable 
during the entire run. The liquid-argon response was monitored with radioactive sources 
of a and (3 particles throughout the run, as were the gains and pedestals of all readout 
channels. Details can be found in Ref. [24|. 



The ECEM calorimeter provides a measurement of energy and position of the electrons 
from the W and Z boson decays. Due to the fine segmentation of the third layer, we can 
measure the position of the shower centroid with a precision of about 1 mm in the azimuthal 
and radial directions. 

We have studied the response of the ECEM calorimeter to electrons in beam tests P,p5|. 
To reconstruct the electron energy we add the signals observed in each EM layer {i = 
1, . . . ,4) and the first FH layer {i = 5) of an array of 5x5 calorimeter towers, centered on 
the most energetic tower, weighted by a layer-dependent sampling weight Si, 



E = SiGi - 6ec 



(6) 
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To determine the sampling weights we minimize 

where the sum runs over all events, ctem is the resolution given in Eq. |^ and p is the beam 
momentum. We obtain A = 3.74 MeV/ADC count, 6ec = -300 MeV, si = 1.47, sa = 1.00, 
S4 = 1.10, and S5 = 1.67. We arbitrarily fix S3 = 1. The value of 6ec depends on the amount 
of uninstrumented material in front of the calorimeter. The parameters Si to S4 weight the 
four EM layers and S5 the first FH layer. Figure ^ shows the fractional deviation of as a 
function of the beam momentum p. Above 20 GeV the non-linearity is less than 0.1%. 
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FIG. 8. The fractional deviation of the reconstructed electron energy from the beam momen- 
tum from beam tests of an ECEM module. 

The fractional energy resolution can be parameterized as a function of electron energy 
using constant, sampling, and noise terms as 



El \ E ) 



+ 



(8) 



with cem = 0.003, sem = 0.157 GeV^/^ and uem = 0.29 GeV in the end calorimeters, as 
measured in beam tests RPS . 



4- Muon Spectrometer 



The D0 muon spectrometer consists of five separate solid-iron toroidal magnets, together 
with sets of proportional drift tube chambers to measure the track coordinates. The central 
toroid covers the region \ri\< 1, two end toroids cover 1 <| |< 2.5, and the small-angle 
muon system covers 2.5 <| |< 3.6. There is one layer of chambers inside the toroids and 
two layers outside for detecting and reconstructing the trajectory and the momentum of 
muons. 



5. Luminosity Monitor 



Two arrays of scintillator hodoscopes, mounted in front of the EC cryostats, register 
hits with a 220 ps time resolution. They serve to detect the occurance of an inelastic pp 
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interaction. The particles from the breakup of the proton give rise to hits in the hodoscopes 
on one side of the detector that are tightly clustered in time. For events with a single 
interaction, the location of the interaction vertex can be determined with a resolution of 
3 cm from the time difference between the hits on the two sides of the detector for use in 
the Level 2 trigger. This array is also called the Level trigger because the detection of an 
inelastic pp interaction is required for most triggers. 

6. Trigger 

Readout of the detector is controlled by a two-level trigger system. Level 1 consists 
of an and-or network that can be programmed to trigger on a pp crossing if a number of 
preselected conditions are satisfied. The Level 1 trigger decision is taken within the 3.5 
/is time interval between crossings. As an extension to Level 1, a trigger processor (Level 
1.5) may be invoked to execute simple algorithms on the limited information available at 
the time of a Level 1 accept. For electrons, the processor uses the energy deposits in each 
trigger tower as inputs. The detector cannot accept any triggers until the Level 1.5 processor 
completes execution and accepts or rejects the event. 

Level 2 of the trigger consists of a farm of 48 VAXstation 4000's. At this level, the 
complete event is available. More sophisticated algorithms refine the trigger decisions and 
events are accepted based on preprogrammed conditions. Events accepted by Level 2 are 
written to magnetic tape for offline reconstruction. 

IV. DATA SELECTION 
A. Trigger 

The conditions required at trigger Level 1 for W and Z boson candidates are: 

• pp interaction: Level hodoscopes register hits consistent with a pp interaction. Using 
monitor trigger data, the efficiency of this condition has been measured to be 98.6%. 

• Main Ring \'cto: No Main Ring proton bunch passes through the detector within 800 
ns of the pp crossing and no protons were injected into the Main Ring less than 400 
ms before the pp crossing. 

• EM trigger towers: There are one or more EM trigger towers with £^ sin ^ > T, where 
E is the energy measured in the tower, 9 is the polar angle of the tower with the beam 
measured from the center of the detector, and T is a programmable threshold. This 
requirement is fully efficient for electrons with pt > 2T. 

The Level 1.5 processor recomputes the transverse electron energy by adding the adjacent 
EM trigger tower with the largest signal to the EM trigger tower that exceeded the Level 1 
threshold. In addition, the signal in the EM trigger tower that exceeded the Level 1 threshold 
must constitute at least 85% of the signal registered in this tower if the hadronic layers are 
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also included. This EM fraction requirement is fully efficient for electron candidates that 



pass our offline selection (Sec. [IV D|) . 

Level 2 uses the EM trigger tower that exceeded the Level 1 threshold as a starting point. 
The Level 2 algorithm finds the most energetic of the four calorimeter towers that make up 
the trigger tower, and sums the energy in the EM sections of a 3x3 array of calorimeter 
towers around it. It checks the longitudinal shower shape by applying cuts on the fraction 
of the energy in the different EM layers. The transverse shower shape is characterized by 
the energy deposition pattern in the third EM layer. The difference between the energies in 
concentric regions covering 0.25x0.25 and 0.15x0.15 in r] x (p must be consistent with an 
electron. Level 2 also imposes an isolation condition requiring 



Y.i Ej sin Oj - pt 
Pt 



< 0.15 



(9) 



where Ei and 6i are the energy and polar angle of cell i, the sum runs over all cells within a 
cone of radius R = a/A0^ + At]'^ = 0.4 around the electron direction and pt is the transverse 
momentum of the electron p6|] . 

The Pt of the electron computed at Level 2 is based on its energy and the 2;-position 
of the interaction vertex measured by the Level hodoscopes. Level 2 accepts events that 
have a minimum number of EM clusters that satisfy the shape cuts and have pt above a 
preprogrammed threshold. Figure |^ shows the measured relative efficiency of the Level 2 
electron filter for forward electrons versus electron p^ for a Level 2 p-r threshold of 20 GeV. 
We determine this efficiency using Z boson data taken with a lower threshold value (16 GeV) 
for one electron. The efficiency is the fraction of electrons above a Level 2 pt threshold of 
20 GeV. The curve is the parameterization used in the fast Monte Carlo (see Sec. 0). 




24 26 28 30 

PT(e) (GeV) 

FIG. 9. The relative efficiency of the Level 2 electron filter for a threshold of 20 GeV for EC 
electrons, as a function of the pt{^) computed offline for the W boson mass analysis. 



Level 2 also computes the missing transverse momentum based on the energy registered 
in each calorimeter cell and the vertex 2;-position as measured by the Level hodoscopes. 
The level 2 W boson trigger requires minimum of 15 GeV. We determine the efficiency 



16 



curve for a 15 GeV Level 2 j^x requirement from data taken without the Level 2 jfj- condition. 
Figure shows the measured efficiency versus Pt(^) as computed for the W mass analysis, 
when the electron is detected in the end calorimeters. The curve is the parameterization 
used in the fast Monte Carlo. 
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FIG. 10. The efficiency of a 15 GeV Level 2 requirement for EC electrons, as a function of 
the pri'^) computed for the W boson mass analysis. 



B. Reconstruction 

1 . Electron 

We identify electrons as clusters of adjacent calorimeter cells with significant energy 
deposits. Only clusters with at least 90% of their energy in the EM section and at least 60% 
of their energy in the most energetic calorimeter tower are considered as electron candidates. 
For most electrons we also reconstruct a track in the CDC or FDC that points towards the 
centroid of the cluster. 

We compute the forward electron energy E{e) from the signals in all cells of the EM 
layers and the first FH layer whose centers lie within a projective cone of radius 20 cm 
and centered at the cluster centroid. In the computation we use the sampling weights 
and calibration constants determined using the test-beam data (Sec. |111B3| ), except for 
the overall energy scale A and the offset ^ec, which we take from an in situ calibration 

(Sec. |vrq ). 

The calorimeter shower centroid position (xcai, Vcai, ^cai), the track coordinates (xtrk, Vtrk, 
Ztrk), and the proton beam trajectory define the electron angle. We determine the position 
of the electron shower centroid Xcai = (a^cai, Z/cai, -2^cai) in the calorimeter from the energy 
depositions in the third EM layer by computing the weighted mean of the positions Xi of 
the cell centers. 
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The weights are given by 



Xcal = ■ (10) 




Wi = max 0, Wo + log ^-Ar- , (11) 



where is the energy in cell i, Wq is a parameter which depends upon ri{e), and E{e) is the 
energy of the electron. The FDC track coordinates are reported at a fixed z position using 
a straight line fit to all the drift chamber hits on the track. The calibration of the radial 
coordinates measured in the cylindrical coordinate system contributes a systematic uncer- 
tainty to the W boson mass measurement. Using tracks from many events reconstructed 
in the vertex drift chamber, we measure the beam trajectory for every run. The closest 
approach to the beam trajectory of the line through the shower centroid and the track coor- 
dinates defines the ^-position of the interaction vertex (^vtx)- The beam trajectory provides 
(3^vtx-?/vtx)- In Z — * ee events, we may have two electron candidates with tracks. In this 
case we take the point determined from the more central electron as the interaction vertex, 
because this gives better resolution. Using only the electron track to determine the position 
of the interaction vertex, rather than all tracks in the event, makes the resolution of this 
measurement less sensitive to the luminosity and avoids confusion between vertices in events 
with more than one pp interaction. 

We then define the azimuth 0(e) and the polar angle 9{e) of the electron using the vertex 
and the shower centroid positions 

tan0(e) = ^^^^^^ , (12) 



tan6'(e) = . (13) 

^ca\ ^vtx 

Neglecting the electron mass, the momentum of the electron is given by 

/ sin 6{e) cos (/ 

p{e)=E{e) sin ^(e) sin 0(e) |. (14) 
\ cos6'(e) 



2. Recoil 

We reconstruct the transverse momentum of all particles recoiling against the W or Z 
boson by taking the vector sum 

^=j:E.sint,.[2^), (15) 

where the sum runs over all calorimeter cells that were read out, except those that belong 
to electron cones. Ei are the cell energies, and (pi and 6i are the azimuth and polar angle of 
the center of cell i with respect to the interaction vertex. 
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3. Derived Quantities 
In the case of Z — ee decays, we define the dielectron momentum 

p{ee) = p(ei) +p(e2) 
and the dielectron invariant mass 



m(eei 



2E{ei)E{e2)il - cosu) 



(16) 



(17) 



where u is the opening angle between the two electrons. It is useful to define a coordinate 
system in the plane transverse to the beam that depends only on the electron directions. We 
follow the conventions first introduced by UA2 [|l^] and call the axis along the inner bisector 
of the transverse directions of the two electrons the ry-axis and the axis perpendicular to 
that the ^-axis. Projections on these axes are denoted with subscripts r] or ^. Figure pi] 
illustrates these definitions. 




FIG. 11. Illustration of momentum vectors in the transverse plane for Z ^ ee candidates. 
The vectors drawn as thick lines are directly measured. 

In the case oiW ^ ev decays, we define the transverse neutrino momentum 

Pj,{i>) = -pT{e) - ut (18) 

and the transverse mass (Eq. |^). Useful quantities are the projection of the transverse recoil 
momentum on the transverse component of the electron direction, 

u\\ = Ut -pTie) , (19) 

and the projection perpendicular to the transverse component of the electron direction, 

u± = ut- {pT{e) X z) . (20) 



Figure O illustrates these definitions. 
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FIG. 12. Illustration of momentum vectors in the transverse plane for W ^ eu candidates. 
The vectors drawn as thick lines are directly measured. 

C. Electron Identification 

1. Fiducial Cuts 

Electrons in the ECEM are defined by the pseudorapidity r] of the cluster centroid 
position with respect to the center of the detector. We define forward electrons by 
1.5 <| r7det(e) |< 2.5. 



2. Quality Variables 

We test how well the shape of a cluster agrees with that expected for an electromagnetic 
shower by computing a quality variable (x^) for all cell energies using a 41-dimensional co- 



variance matrix. The covariance matrix was determined from GEANT-based simulations 
that were tuned to agree with extensive test beam measurements. 



To determine how well a track matches a cluster, we extrapolate the track to the third 
EM layer in the end calorimeter and compute the distance between the extrapolated track 
and the cluster centroid in the azimuthal direction, As, and in the radial direction, Ap. The 
variable 




<k= TT + ^ (21) 



quantifies the quality of the match. The parameters 5s = 0.25 cm and 6p = 1.0 cm are the 
resolutions with which As and Ap are measured, as determined using the end calorimeter 
electrons from W ^ ev decays. 

In the EC, electrons must have a matched track in the forward drift chamber to suppress 
background due to misidentification. In the CC, we define "tight" and "loose" criteria. The 
tight criteria require a matched track in the CDC, defined as the track with the smallest 
(Jtrk- The loose criteria do not require a matched track and help increase the electron finding 
efficiency for Z — > ee decays with at least one central electron. 
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The isolation fraction is defined as 



Jiso = ^ , [^^) 



^core 



where i?cone is the energy in a cone of radius R = 0.4 around the direction of the electron, 
summed over the entire depth of the calorimeter, and -Ecore is the energy in a cone of R = 0.2, 
summed over the EM calorimeter only. 

We use the dE/dx information provided by the FDC on the tracks associated with the 
EM calorimeter cluster. The dE/dx information helps to distinguish between singly- ionizing 
electron tracks and doubly-ionizing tracks from photon conversions. 

We identify electron candidates in the forward detectors by making loose cuts on the 
shower shape x^i the track-cluster match quality, and the shower electromagnetic energy 
fraction. The electromagnetic energy fraction is the ratio of the cluster energy measured in 
the electromagnetic calorimeter to the total cluster energy (including the hadronic calorime- 
ter), and is a measure of the longitudinal shower profile. We then use a cut on a 4- variable 
likelihood ratio A4 which combines the information in these variables and the track dE/dx 
into a single variable. The final cut on the likelihood ratio A4 gives the maximum discrimi- 
nation between electrons and jet background, i.e. gives the maximum background rejection 
for any given electron selection efficiency. 



Figure |T3| shows the distributions of the quality variables for electrons in the EC data; 



the arrows indicate the cut values. Table | summarizes the electron selection criteria. 
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FIG. 13. Distributions of the EC electron identification variables for W ^ ev candidates in 
the data. The arrows indicate the cut values. 



D. Data Samples 



The data were collected during the 1994-1995 Tevatron run. After the removal of runs 
in which parts of the detector were not operating adequately, the data correspond to an 
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TABLE I. Electron selection criteria. A^^ai is the difference in azimuthal angle between the 
cluster centroid and the CC module edge. 



variable 



CC (loose) 
|A(/)eal| > 0.02 
I Zeal I < 108 cm 



< 100 
/iso < 0.15 



CC (tight) 



|A(/)cal| > 0.02 
I Zeal I < 108 cm 
l^trkl < 80 cm 
< 100 
/iso < 0.15 

fJtrk < 5 



EC (tight) 



fiducial cuts 



shower shape 
isolation 
track match 
4- variable 
likelihood ratio 



1.5 < \r]\ < 2.5 

< 200 
/i.o < 0.15 

<7trk < 10 

A4 <4 



integrated luminosity of 82 pb^^ We select W boson decay candidates by requiring: 
Level 1: pp interaction 

Main Ring Veto 

EM trigger tower above 10 GeV 
Level 1.5: > 1 EM cluster above 15 GeV 
Level 2: electron candidate with pr > 20 GeV 

momentum imbalance > 15 GeV 
offline: > 1 tight electron candidate in EC 

PT{e) > 30 GeV 

PtIi^) > 30 GeV 

ut < 15 GeV 

This selection gives us 11,089 W boson candidates. We select Z boson decay candidates by 
requiring: 

Level 1: pp interaction 

> 2 EM trigger towers above 7 GeV 

Level 1.5: > 1 EM cluster above 10 GeV 

Level 2: > 2 electron candidates with p^ > 20 GeV 

offline: > 2 electron candidates 
pr(e) > 30 GeV (EC) 
or pT{e) > 25 GeV (CC) 
We accept Z ^ ee decays with at least one electron candidate in the EC and the other in the 
CC or the EC. EC candidates must pass the tight electron selection criteria. A CC candidate 
may pass only the loose criteria. We use the 1,687 events with at least one electron in the EC 
(CC/EC + EC/EC Z samples) to calibrate the calorimeter response to electrons (Sec. |V|). 
These events need not pass the Main Ring Veto cut because Main Ring background does 
not affect the EM calorimeter. Of these events, those that do pass the Main Ring Veto have 
been used to calibrate the recoil momentum response. The events for which both electrons 
are in the EC (EC/EC Z sample) and which pass the Main Ring Veto serve to check the 
calibration of the recoil response (Sec. |V11| ). Table |1| summarizes the data samples. 
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TABLE II. Number of W and Z boson candidate events. 



channel 




Z ^ ee 




W ^eu 


fiducial region of electrons 


CC/EC 




EC/EC 


EC 




1265 




422 


11089 



Figure |1J shows the luminosity of the colliding beams during the W and Z boson data 
collection. 
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FIG. 14. The instantaneous luminosity distribution of the W (top) and the Z (bottom) boson 
samples. 

On several occasions we use a sample of 295,000 random pp interaction events for cal- 
ibration purposes. We collected these data concurrently with the W and Z signal data, 
requiring only a pp interaction at Level 1. We refer to these data as "minimum bias events." 

V. FAST MONTE CARLO MODEL 
A. Overview 



The fast Monte Carlo model consists of three parts. First we simulate the production 
of the W OI Z boson by generating the boson four-momentum and other characteristics 
of the event such as the z-position of the interaction vertex and the luminosity. The event 
luminosity is required for luminosity-dependent parametrizations in the detector simulation. 
Then we simulate the decay of the boson. At this point we know the true pt of the boson 
and the momenta of its decay products. We next apply a parameterized detector model 
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to these momenta to simulate the observed transverse recoil momentum and the observed 
electron momenta. 

Our fast Monte Carlo program is very similar to the one used in our published CC 
analysis @] , with some modifications in the simulation of forward electron events. 



B. Vector Boson Production 



To specify the production dynamics of vector bosons in pp collisions completely, we 
need to know the differential production cross section in mass Q, rapidity y, and transverse 
momentum qt of the produced W bosons. To speed up the event generation, we factorize 
this into 



a 



d^a 



dqj^dydQ dqj^dy 



X 



-Ml 



da 
dQ 



(23) 



to generate qr, y, and Q of the bosons. 

For pp collisions, the vector boson production cross section is given by the parton cross 
section dij convoluted with the parton distribution functions (pdf) f{x,Q^) and summed 
over parton flavors 



d'a 
dq^dy 



j dxi j dx2fi{xi,Q'^)fj{x2,Q'^) 

id 



2\ d aij 



6{sXiX2 - Q J , 2 A 

dqt^dy 



(24) 



The cross section d a / dqj^dy\Q2=M'^ has been computed by several authors p9| , |30| using a 



perturbative calculation |^ for the high-g-p regime and the Collins-Soper resummation for- 
malism [0,0 for the low-gr regime. We use the code provided by the authors of Ref. ||2^ 
and the MRST parton distribution functions p3 to compute the cross section. The produc- 
tion of WW , WZ and is suppressed by three orders of magnitude compared to inclusive 
W production. 

We use a Breit-Wigner curve with a mass-dependent width for the line shape of the W 
boson. The intrinsic width of the W is Ty/ = 2.062 ± 0.059 GeV The line shape is 
skewed due to the momentum distribution of the quarks inside the proton and antiproton. 
The mass spectrum is given by 



da 
dQ 



,2 , Q'^w 



(25) 



We call 



A.(Q) = — E 



1 dx 

)2/s X 



Mx,Q')f,{Qysx,Q' 



(26) 



the parton luminosity. To evaluate it, we generate W ^ ev events using the HERWIG Monte 
Carlo event generator pH], interfaced with pdflib p^, and select the events subject to the 



24 



same fiducial cuts as for tlie W and Z boson samples with at least one electron in EC. We 
plot the mass spectrum divided by the intrinsic line shape of the W boson. The result is 
proportional to the parton luminosity, and we parameterize the shape of the spectrum with 
the function ^ 

^M = -Q- ■ (27) 

Table ^ shows the parton luminosity slope /3 for W and Z events for the different topologies. 
The value of (3 depends on the rapidity distribution of the W and Z bosons, which is 
restricted by the fiducial cuts that we impose on the decay leptons. The values of (3 given in 
Table |IT| are for the rapidity distributions of W and Z bosons that satisfy the fiducial cuts 



given in Sec. I^. The uncertainty in P is about 0.001 GeV , due to Monte Carlo statistics 



and uncertainties in the acceptance. 



TABLE III. Parton luminosity slope f3 in the W and Z boson production model. The /? value 
is given for W ^ ev decays with the electron in the EC and for Z ^ ee decays with at least one 
electron in the EC. 





Z production 


W production 




/3 (GeV-i) 


/3 (GeV-i) 


CC/EC 


9.9 X 10-3 




EC/EC 


19.9 X 10-^ 




EC 




16.9 X 10-3 



Bosons can be produced by the annihilation of two valence quarks, two sea quarks, or 
one valence quark and one sea quark. Using the HERWIG events, we evaluate the fraction /gg 
of bosons produced by the annihilation of two sea quarks. We find fss = 0.207, independent 
of the boson topology. 

To generate the boson four-momenta, we treat da/dQ and d'^a/dq^dy as probability 
density functions and pick Q from the former and a pair of y and qx values from the latter. 
For a fraction fss the boson helicity is +1 or —1 with equal probability. The remaining W 
bosons always have helicity —1. Finally, we pick the z-position of the interaction vertex 
from a Gaussian distribution centered at z = with a standard deviation of 27 cm and a 



luminosity for each event from the histogram in Fig. 14 



C. Vector Boson Decay 

At lowest order, the boson is fully polarized along the beam direction due to the 
V ^ A coupling of the charged current. The resulting angular distribution of the charged 
lepton in the W boson rest frame is given by 

'^'^ oc (l-Agcosr)2 , (28) 



d cos 9* 
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where A is the hehcity of the W boson with respect to the proton direction, q is the charge 
of the lepton, and 6* is the angle between the charged lepton and proton beam directions 
in the W rest frame. The spin of the W boson points along the direction of the incoming 
antiquark. Most of the time, the quark comes from the proton and the antiquark from the 
antiproton, so that A = — 1. Only if both quark and antiquark come from the sea of the 
proton and antiproton, is there a 50% chance that the quark comes from the antiproton and 
the antiquark from the proton and in that case A = 1 (see Fig. ITH). 



P P P P 

FIG. 15. Polarization of the W boson produced in pp collisions if the quark comes from the 
proton (left) and if the antiquark comes from the proton (right). The short thick arrows indicate 
the orientations of the particle spins. 

When 0{as) processes are included, the boson acquires finite transverse momentum and 
Eq. H becomes m 

oc (l - Xqaiiqr) cos 6*03 + a2(gr) cos^ 6*08) (29) 



d cos 6*08 

for W bosons after integration over 0. The angle ^C8 in Eq. is now defined in the Collins 



Soper frame The values of ai and 02 as a function of transverse boson momentum 

have been calculated at C(«s) [0. We have implemented the angular distribution given in 
Eq. ^ in the fast Monte Carlo. The angular distribution of the leptons from Z ^ ee decays 
is also generated according to Eq. ^ but with ai and 02 computed for Z ^ ee decays . 



Radiation from the decay electron or the W boson biases the mass measurement. If the 
decay electron radiates a photon and the photon is sufficiently separated from the electron 
so that its energy is not included in the electron energy, or if an on-shell W boson radiates 
a photon and therefore is off-shell when it decays, the measured mass is biased low. We use 



the calculation of Ref. |^ to generate W ewj and Z eej decays. The calculation 
gives the fraction of events in which a photon with energy -^(7) > Eq is radiated, and the 
angular distribution and energy spectrum of the photons. Only radiation from the decay 
electron and the W boson, if the final state W is off-shell, is included to order a. Radiation 
by the initial quarks or the W boson, if the final W is on-shell, does not affect the mass 
of the eu pair from the W decay. We use a minimum photon energy Eq = 50 MeV, and 
calculate that in 30.6% of all W decays a photon with -£^(7) > 50 MeV is radiated. Most of 
these photons are emitted close to the electron direction and cannot be separated from the 
electron in the calorimeter. For Z — »• ee decays, there is a 66% probability for either of the 
electrons to radiate a photon with £^(7) > 50 MeV. 

If the photon and electron are close together, they cannot be separated in the calorime- 
ter. The momentum of a photon with AR{e'~f) < Rq is therefore added to the electron 
momentum, while for AR{e'y) > Rq, a photon is considered separated from the electron and 
its momentum is added to the recoil momentum. We use Rq = 20 cm, which is the size of 
the cone in which the electron energy is measured. We refer to i?o as the photon coalescing 
radius. 
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W boson decays through the channel W ^ tv ^ evvv are topologically indistinguish- 
able from W ^ ev decays. We therefore include these decays in the W decay model, properly 
accounting for the polarization of the tau leptons in the decay angular distributions. In the 
standard model and neglecting small phase space effects, the fraction of W boson decays to 
electrons that proceed via tau decay is B{t euV)/ (1 + B{t euV)) = 0.151. 



D. Detector Model 

The detector simulation uses a parameterized model for detector response and resolution 
to obtain a prediction for the distributions of the observed electron and recoil momenta. 

When simulating the detector response to an electron of energy Eq, we compute the 
observed electron energy as 

E{e) = aEcEo + AE{C,r],uii) + aEMX , (30) 

where a^c is the response of the end electromagnetic calorimeter, AE is the energy due to 
particles from the underlying event within the electron cone (parameterized as a function 
of luminosity C, r] and u\\), (Tem is the energy resolution of the electromagnetic calorimeter, 
and X is a random variable from a normal parent distribution with zero mean and unit 
width. 

The transverse energy measurement depends on the measurement of the electron direc- 
tion as well. We determine the shower centroid position by intersecting the line defined 
by the event vertex and the electron direction with a plane perpendicular to the beam and 
located at z = ± 179 cm (the longitudinal center of the ECEM3 layer). We then smear 
the azimuthal and radial coordinates of the intersection point by their resolutions. We de- 
termine the radial coordinate of the FDC track by intersecting the same line with a plane 
at z = ±105 cm, the defined z position of the FDC track centroid, and smearing by the 
resolution. The measured angles are then obtained from the smeared points as described in 



Section IV Bl 



The model for the particles recoiling against the W boson has two components: a "hard" 
component that models the pt of the W boson, and a "soft" component that models detector 
noise and pile-up. Pile-up refers to the effects of additional pp interactions in the same or 
previous beam crossings. For the soft component we use the transverse momentum balance 

measured in minimum bias events recorded in the detector. The minimum bias events 
are weighted so that their luminosity distribution is the same as that of the W sample. The 
observed recoil pr is then given by 

-Au\\{C,r],u\()pT{e) 

+Oimb§T , (31) 

where qt is the generated value of the boson transverse momentum, Rrec is the (in general 
momentum-dependent) response, cJrec is the resolution of the calorimeter (parameterized as 
c"rec = -SrecA/wr), Au\\ is the trausvcrsc energy flow into the electron window (parameterized 
as a function of £, r) and uii), and a^h is a correction factor that allows us to adjust the 
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resolution to the data, accounting for the difference between the data minimum bias events 
and the underlying spectator collisions in W events. The quantity Au\\ is different from the 
transverse energy added to the electron, AEt, because of the difference in the algorithms 
used to compute the electron Et and the recoil pr- 

We simulate selection biases due to the trigger requirements and the electron isolation 
by accepting events with the estimated efficiencies. Finally, we compute all the derived 
quantities from these observables and apply fiducial and kinematic cuts. 



VI. ELECTRON MEASUREMENT 



A. Angular Calibrations 



The FDC detectors have been studied and calibrated extensively in a test beam [^]. We 



use collider data muons which traverse the forward muon detectors and the FDC to provide 
a cross-check of the test beam calibration of the radial measurement of the track in the 
FDC. We predict the trajectory of the muon through the FDC by connecting the hits in the 
innermost muon chambers with the reconstructed event vertex by a straight line. The FDC 



track coordinate can then be compared relative to this line. Figure |I6| shows the difference 
between the predicted and the actual radial positions of the track. These data are fit to a 
straight line constrained to pass through the origin. We find the track position is consistent 
with the predicted position. 
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FIG. 16. Residue of the radial position of the FDC track centroid from the predicted radial 
position of forward muon tracks at the FDC, as a function of the track radial position. The solid 
line is a fitted straight line constrained to pass through the origin. 



We calibrate the shower centroid algorithm using Monte Carlo electrons simulated using 
GEANT and electrons from the Z ^ ee data. We apply a polynomial correction as a function 
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of Teal and the distance from the cell edges based on the Monte Carlo electrons. We refine the 
calibration with the Z ^ ee data by exploiting the fact that both electrons originate from 



the same vertex. Using the algorithm described in Sec. B 1| , we determine a vertex for each 
electron from the shower centroid and the track coordinates. We minimize the difference 



between the two vertex positions as a function of an rcai scale factor /5ec (see Fig. [T^). The 
correction factor is /3ec = 0.9997 ± 0.00044 for EC North, and /3ec = 1-00225 ± 0.00044 for 
EC South. We find no systematic radial dependence of these correction factors. 

We quantify the FDC and EC radial calibration uncertainty in terms of scale factor 
uncertainties 5/3fdc = ±0.00054 and (5/5ec = ±0.0003 for the radial coordinate. The un- 
certainties in these scale factors lead to a 20 MeV uncertainty in the EC W boson mass 
measurement. 



581r 



580' 

0.999 



586r 



585' 

1.001 1.001 



Pgp(north) 



T 



1.003 

PgpCsouth) 



FIG. 17. The versus /?ec value. 



B. Angular Resolutions 

The resolution for the radial coordinate of the track, rtrk, is determined from the Z 
ee sample. Both electrons originate from the same interaction vertex and therefore the 
difference between the interaction vertices reconstructed from the two electrons separately, 
^vtx(ei) — -Zvtx(e2), is a measure of the resolution with which the electrons point back to the 
vertex. The points in Fig. ^ show the distribution of z^t^{ei) — -2vtx(e2) observed in the 
CC/EC and EC/EC Z samples with matching tracks required for both electrons. 
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FIG. 18 

Z ^ ee samples (•) and the fast Monte Carlo simulation ( 



A Monte Carlo study based on single electrons generated with a geant simulation shows 
that the resolution of the shower centroid algorithm is 0.1 cm in the EC, consistent with EC 
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electron beam tests. We then tune the resolution function for rtrk in the fast Monte Carlo 
so that it reproduces the shape of the Zvtx(ei) — Zvtx{e2) distribution observed in the data. 
We find that a resolution function consisting of two Gaussians 0.2 cm and 1.7 cm wide, with 
20% of the area under the wider Gaussian, fits the data well. The histogram in Fig. |T8 
shows the Monte Carlo prediction for the best fit, normalized to the same number of events 
as the data. 



C. Underlying Event Energy 

We define a cone which is projective from the center of the detector, has a radius of 
20 cm at the z position of ECEM3 and is centered on the electron cluster centroid. The 
cone extends over the four ECEM layers and the first ECFH layer. This cone contains the 
entire energy deposited by the electron shower plus some energy from other particles. The 
energy in the window is excluded from the computation of ut- This causes a bias in u\i, the 
component of ut along the direction of the electron. We call this bias Au\\. It is equal to 
the momentum fiow observed in the EM and first FH sections of a projective cone of radius 
20 cm at ECEM3. 

We use the W data sample to measure Au\\. For every electron in the W sample, we 
compute the energy fiow into an azimuthally rotated position, keeping the cone radius and 
the radial position the same. For the rotated position we compute the measured transverse 
energy. Since the 7](j) area of the cone increases as the electron rj increases, it is convenient 
to parameterize the transverse energy density, Au\\/6ri6(j). 

At higher luminosity the average number of interactions per event increases and there- 
fore Auw/S'qScj) increases (Fig. |l^). The mean value of Au\\/5ri5(j) increases by 40 MeV 
per 10^°cm~^s~^. The underlying event energy fiow into the electron cone depends on the 
electron 77, as shown in Fig. corrected back to zero luminosity. 

The underlying event energy fiow into the electron cone also depends on the overlap 
between the recoil and the electron. We have found that the best measure of the recoil 
overlap is the component of the total recoil in the direction of the electron, which is u\\. 
Figure |21| shows {Au\\/6ri6(f){C = 0,\r]\= 2.0)), the mean value for Au\\/5r]5(l) corrected to 
zero luminosity and | ?7 1= 2.0, as a function of u\\. In the fast Monte Carlo model, a value 
Au\\/5r]5(l) is picked from the distribution shown in Fig. ^for every event, corrected for u\\, 
rj, and luminosity dependences, and then scaled by the 6ri6(f) area of a 20 cm cone at the 
electron rj. 

The measured electron transverse energy is biased upwards by the additional energy AEt 
in the window from the underlying event. AEt is not equal to Au\\ because the electron 
Et is calculated by scaling the sum of the cell energies by the electron angle, whereas ut is 
obtained by summing the Et of each cell. The ratio of the two corrections as a function of 
electron rj is shown in Fig. |23|. 

The uncertainty in the underlying event transverse energy density has a statistical com- 
ponent (14 MeV) and a systematic component (24 MeV). The systematic component is 
derived from the difference between the measurement close to the electron (where it is bi- 
ased by the isolation requirement) and far from the electron (where it is not biased). The 
total uncertainty in the underlying event transverse energy density is 28 MeV. 
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FIG. 19. The instantaneous luminosity dependence of {l^u\\/5rj5(j)). 




lTl(e)l 

FIG. 20. The variation of {Auii/5ri6(j)) as a function of electron rj. 
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FIG. 21. The variation of {Au\\/ 5r]5(j)) as a function of uy. The region between the arrows is 
populated by the W boson sample. 
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FIG. 22. The distribution of /\u\\/ drjdcj) in the W signal sample, corrected to £=0, | r] 
wii=0 . 
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FIG. 23. The ratio of the {l!^u\\/ 5r]6(j)) corrections to the electron and the recoil as a function 
of electron rj. 



D. U|| Efficiency 

The efficiency for electron identification depends on the electron environment. Well- 
isolated electrons are identified correctly more often than electrons near other particles. 
Therefore W decays in which the electron is emitted in the same direction as the particles 
recoiling against the W boson are selected less often than W decays in which the electron 
is emitted in the direction opposite the recoiling particles. This causes a bias in the lepton 
Pt distributions, shifting Pt{g) to larger values and Pt{i^) to lower values, whereas the rriT 
distribution is only shghtly affected. 

We measure the electron finding efficiency as a function of u\\ using Z ^ ee events. The Z 
event is tagged with one electron, and the other electron provides an unbiased measurement 
of the efficiency. Following background subtraction, the measured efficiency is shown in 
Fig. ^ The line is a fit to a function of the form 



Bo 



S{U\\ - Mo) 



for u\\ < Uq 
otherwise. 



(32) 



The parameter Eq is an overall efficiency which is inconsequential for the W mass measure- 
ment, Uq is the value of at which the efficiency starts to decrease as a function of u\\, and 
s is the rate of decrease. We obtain the best fit for uq = —2.4 GeV and s = 0.0029 GeV~^. 
These two values are strongly anti-correlated. The error on the slope 6s = ±0.0012 GeV~^ 
accounts for the statistics of the Z sample. 
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FIG. 24. The EC electron selection efficiency as a function of uy. 

E. Electron Energy Response 

Equation |^ relates the reconstructed electron energy to the recorded end calorimeter 
signals. Since the values for the constants were determined in the test beam, we determine 
the offset Sec and a scale oec, which essentially modifies A, in situ with collider Z ^ ee 
data. 

The electrons from Z decays are not monoenergetic and therefore we can make use of 
their energy spread to constrain 6ec- When both electrons are in the EC, we can write 

m(ee) = oecMz + IzSec (33) 

for Sec ^ E{ei) + £"(62). fz is a kinematic function related to the boost of the Z boson, 
and is given by fz = [E{ei) +E{e2)]{l — cosu)/m{ee), where u is the opening angle between 
the two electrons. When one electron is in the CC and one is in the EC, we can write 

m{ee) = a/«cc«ecA^z + /z^ec, (34) 

where fz = E{e2){l — cos u) / m{ee) and 62 is the CC electron. When we apply this formula, 
we have already corrected the CC electron for the corresponding CCEM offset, ^cc = —0.16 
GeV, which was measured for our CC W mass analysis acc is the CC electromagnetic 
energy scale, which is determined by fitting the m(ee) spectrum of the CC/CC Z sample. 

We plot m(ee) versus fz and extract Sec as the slope of the fitted straight line. We use 
the fast Monte Carlo to correct for residual biases introduced by the kinematic cuts. The 
Sec measurements from the CC/EC and EC/EC Z samples are shown in Fig. along with 
the statistical uncertainties. We obtain the average Sec = —0.1 ±0.7 GeV. The uncertainty 
in this measurement of ^ec is dominated by the statistical uncertainty due to the finite 
size of the Z sample. As Fig. ^ shows, the offsets measured in the north and south end 
calorimeters separately are completely consistent. 
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FIG. 25. The ECEM offset measurements using the CC/EC and EC/EC Z samples. The 
labels indicate the calorimeter cryostat in which each of the Z decay electrons was detected. CC 
indicated the central calorimeter and ECN (ECS) indicates the north (south) end calorimeter 
respectively. 



After correcting the data with this value of 5ec we determine Q!ec so that the position of 
the Z peak predicted by the fast Monte Carlo agrees with the data. To determine the scale 
factor that best fits the data, we perform a maximum likelihood fit to the m(ee) spectrum 
between 70 GeV and 110 GeV. In the resolution function we allow for background shapes 
determined from samples of events with two EM clusters that fail the electron quality cuts 
(Fig. |23). The background normalization is obtained from the sidebands of the Z peak. 
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FIG. 26. The dielectron mass spectrum from the CC/EC (left) and EC/EC (right) samples 
of events with two EM clusters that fail the electron quality cuts. The superimposed curves shows 
the fitted functions used to model the shape of the background in the Z samples. 

Figure ^ shows the m(ee) spectrum for the CC/EC Z sample and the Monte Carlo 
spectrum that best fits the data for (5ec = —0.1 GeV. The for the best fit to the CC/EC 
m(ee) spectrum is 14 for 19 degrees of freedom. For age = 0.95143 ± 0.00259, the Z peak 
position of the CC/EC sample is consistent with the known Z boson mass. The error reflects 
the statistical uncertainty. The background has no measurable effect on the result. 
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FIG. 27. The dielectron mass spectrum from the CC/EC Z sample. The superimposed curve 
shows the maximum hkehhood fit and the shaded region the fitted background. 



Figure ^ shows the m(ee) spectrum for the EC/EC Z sample and the Monte Carlo 
spectrum that best fits the data for 5ec = —0.1 GeV. The for the best fit to the EC/EC 
m(ee) spectrum is 12 for 17 degrees of freedom. For oec = 0.95230 ± 0.00231, the Z peak 
position of the EC/EC sample is consistent with the known Z boson mass. The error reflects 
the statistical uncertainty and the uncertainty in the background. 
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FIG. 28. The dielectron mass spectrum from the EG/EC Z sample. The superimposed curve 
shows the maximum likelihood fit and the shaded region the fitted background. 

Combining the oec measurements from the CC/EC and the EC/EC Z samples, we 
obtain the ECEM energy scale 



ftEC = 0.95179 ±0.00187 



(35) 
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The difference between the ECEM scales measured separately in the north and south 
calorimeters is 0.0040 ± 0.0037, consistent with the calorimeters having the same EM re- 
sponse. 

F. Electron Energy Resolution 

Equation § gives the functional form of the electron energy resolution. We take the 
intrinsic resolution of the end calorimeter, which is given by the sampling term sem, from 
the test beam measurements. The noise term uem is represented by the width of the electron 
underlying event energy distribution (Fig. 0). We measure the constant term cem from the 
Z line shape of the data. We fit a Breit-Wigner convoluted with a Gaussian, whose width 
characterizes the dielectron mass resolution, to the Z peaks for the CC/EC and EC/EC 
samples separately. Figure ^ shows the width am(ee) of the Gaussian fitted to the Z peak 
predicted by the fast Monte Carlo as a function of cem- The horizontal lines indicate 
the width of the Gaussian fitted to the Z samples and its uncertainties. For the data 
measurements of 

am = 2.47 ± 0.05 GeV (CC/EC) 

am = 2.72 ± 0.11 GeV (EC/EC) (36) 

we extract from the CC/EC Z boson events cec = ^-Q-i's^o and from the EC/EC Z events 
we extract cec = O.OtQ;Q%. We take the combined measurement to be 

CEC = i.ot?:^ %■ (37) 

The measured Z boson mass does not depend on cec- 

VII. RECOIL MEASUREMENT 
A. Recoil Momentum Response 

The detector response and resolution for particles recoiling against a W boson should be 
the same as for particles recoiling against a Z boson. For Z — > ee events, we can measure 
the transverse momentum of the Z boson from the e+e~ pair, ^^(ee), into which it decays, 
and from the recoil momentum ut in the same way as for W ^ eu events. By comparing 
Pri^e) and ut, we calibrate the recoil response relative to the electron response. 

The recoil momentum is carried by many particles, mostly hadrons, with a wide momen- 
tum spectrum. Since the response of the calorimeter to hadrons is slightly nonlinear at low 
energies, and the recoil particles see a reduced response at module boundaries, we expect 
a momentum- dependent response function with values below unity. To fix the functional 
form of the recoil momentum response, we studied the response predicted by a Monte 
Carlo Z ^ ee sample obtained using the herwig program and a GEANT-based detector sim- 
ulation. We projected the reconstructed transverse recoil momentum onto the transverse 
direction of motion of the Z boson and define the response as 
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FIG. 29. The dielectron mass resolution versus the constant term cem- The top plot is for the 
CC/EC Z events and the bottom plot is for the EC/EC Z events. 



(38) 



where qr is the generated transverse momentum of the Z boson. A response function of the 
form 



Rrec = ttrec + /^rec In [flT j GeV) 



(39) 



fits the response predicted by geant with Orec = 0.713 ± 0.006 and /?rec = 0.046 ± 0.002. 



This functional form also describes the jet energy response |^ of the D0 calorimeter. 

The recoil response for data was calibrated against the electron response by requiring j)^ 
balance in Z — > ee decays for our published CC analysis . The Z boson measured with 
the electrons and the recoil are projected on the r/ axis, defined as the bisector of the two 
electron directions in the transverse plane. From the CC/CC + CC/EC Z boson events, 
we measured ttrec = 0.693 ± 0.060 and /^rec = 0.040 ± 0.021, in good agreement with the 
Monte Carlo prediction. To compare the recoil response measured with Z events of different 
topologies, we scale the recoil measurement with the inverse of the response parametrization 



i?,ec = 0.693 + 0.04 ■ In (pT(ee)/GeV) 



(40) 



and plot the sum of the projections versus p^(ee), as shown in Fig. We see no Pr^iee) 
dependence to the balance measured using the Z boson events with at least one central 
electron, since this sample was used to derive the values of these parameters. The EC/EC 
Z boson events give a recoil response measurement statistically consistent with the above. 
Hence we use the same recoil response for the EC and the CC W boson events H]. 
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FIG. 30. The recoil momentum response in the CC/CC + CC/EC (left) and the EC/EC 
(right) Z samples as a function of pr^{ee). 



B. Recoil Momentum Resolution 



The widths of the balance and the balance (where the ^ axis is perpendicular to the 
7] axis) are sensitive to the recoil resolution. Figures |3ll-p^ show the comparison between 
the data and Monte Carlo for the recoil resolution determined in our CC W mass analysis 
1^. The balance width is in good agreement between data and Monte Carlo for all Z 
boson topologies. Hence we use the same recoil resolution for EC W boson events as for the 
CC W boson events [g. 
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FIG. 31. The r/-balance distribution for the Z boson data (•) and the fast Monte Carlo 

simulation ( ). The plot on the left is for the CC/CC + CC/EC Z events and the plot on the 

right is for the EC/EC Z events. 
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FIG. 32. The ^-balance distribution for the Z boson data (•) and the fast Monte Carlo 

simulation ( ). The plot on the left is for the CC/CC + CC/EC Z events and the plot on the 

right is for the EC/EC Z events. 
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C. Comparison with W Boson Data 



We compare the recoil momentum distributions in the W boson data to the predictions 
of the fast Monte Carlo, which includes the parameters described in this section and Sec. [V^ . 



Figure ^ shows the u\\ spectra from Monte Carlo and W data. The agreement means that 
the recoil momentum response and resolution and the u\\ efficiency parameterization describe 



the data well. Figures |34|-pq show u_i_, ut, and the azimuthal difference between electron 
and recoil directions from Monte Carlo and W boson data. The figures also show the mean 
and r.m.s. of the data and Monte Carlo distributions and the over the number of degrees 
of freedom (dof). 
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FIG. 33. The spectrum for the W data (•) and the Monte Carlo simulation (- 
mean (/i) and r.m.s. (a) of the distributions and the x^/dof is also shown. 
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VIII. CONSTRAINTS ON THE W BOSON RAPIDITY SPECTRUM 



In principle, if the acceptance for the W —>■ eu decays were complete, the transverse 
mass distribution or the lepton pt distributions would be independent of the W rapidity. 
However, cuts on the electron angle in the laboratory frame cause the observed distributions 
of the transverse momenta to depend on the W rapidity. Hence a constraint on the W 
rapidity distribution is useful in constraining the production model uncertainty on the W 
mass. 

The pseudorapidity distribution of the electron from W ^ eu decays is correlated with 
the rapidity distribution of the W boson. Therefore we can compare the electron t] distri- 
bution between data and Monte Carlo. 

To compare the data with the Monte Carlo, we need to correct for the jet background 
in the data and the electron identification efficiency as a function of rj. We obtain the jet 
background fraction as a function of rj by counting the number of W events that fail electron 
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FIG. 34. The u± spectrum for the W data (•) and the Monte Carlo simulation ( ). 

mean (/i) and r.m.s. (a) of the distributions and the x^/dof is also shown. 
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FIG. 35. The recoil momentum {ut) spectrum for the W data (•) and the Monte Carlo 
simulation ( ). The mean (/i) and r.m.s. (a) of the distributions and the x^/dof is also shown. 
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FIG. 36. The azimuthal difference between electron and recoil directions for the W data (•) 

and the Monte Carlo simulation ( ). The mean (/x) and r.m.s. (a) of the distributions and the 

X^/dof is also shown. 



cuts (see Sec. [iXB| ) in bins of t], subtracting the small contamination due to true electrons, 
and normalizing the entire distribution to the total background fraction (separately in the 
CC and EC). The normalized background r] distribution is subtracted from the t] distribution 
of the data. 

The electron identification efficiency (after fiducial and kinematic cuts) is measured using 
the CC/CC and CC/EC Z ee events. All the electron identification cuts are used to 
identify one electron to tag the event. Candidates are selected in the mass range 81 < rriee < 
101 GeV. Sidebands in the mass range 60 < rriee < 70 GeV and 110 < niee < 120 GeV are 
used for background subtraction. The number of events in which the second electron also 
satisfies all the electron identification cuts is used to calculate the efficiency. The efficiency 



measured in bins of the r] of the second electron is shown in Fig. 37. 
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FIG. 37. Dependence of electron identification efficiency on electron pseudorapidity. Statistical 
errors are shown. 



We scale the electron t] distribution predicted by the Monte Carlo by the r^-dependent effi- 
ciency, and compare to the background-subtracted data in Fig. The errors on the Monte 
Carlo points include the statistical errors on the Monte Carlo sample and the statistical 
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errors on the efficiency measurements. The errors on the data points include the statisti- 
cal errors on the number of candidate events and the statistical errors on the background 



estimate which has been subtracted. Figure 39 shows the ratio between the background- 



subtracted data and the efficiency-corrected Monte Carlo, with the uncertainties mentioned 
above added in quadrature. The Monte Carlo has been normalized to the data. The x^/dof 
shown is with respect to unity. There is good agreement between the data and the Monte 
Carlo. 
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FIG. 38. r] distribution of the electron from TF — > ez/ decays from background-subtracted 
data (•), efficiency-corrected Monte Carlo (o) and the jet background (shaded histogram). The 
distributions drop near |?7|= 1.2 because there is no EM calorimetry in the range 1.1 <|7?det |< 1-4. 




FIG. 39. The ratio of the background-subtracted data and efficiency-corrected Monte Carlo. 
The Monte Carlo has been normalized to the data. The x^/dof is with respect to unity. 

To extract a constraint on the y distribution of the W boson, we introduce in the Monte 
Carlo a scale factor as follows: 

yw ^ krj- yw (41) 

i.e. the rapidity of the W is scaled by the factor /c^. We then compute the between the 
data and Monte Carlo ri{e) distributions for different k^. The result is shown in Fig. HO for 
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the MRS (A') parton distribution functions. Table 
the is minimized for the different pdf 's. 
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FIG. 40. of the electron rj distribution ratio between data and Monte Carlo from unity, as 
a function of the W rapidity scale factor k^. There are 11 degrees of freedom. The Monte Carlo 
uses the MRS (A') parton distribution functions. The horizontal lines indicate Xmm Xmin + 

The uncertainty in is 1.6%, which is the change in A;^ that causes the fo rise by 
one unit above the minimum. We generate Monte Carlo events with different values of /c^ 
and fit them with templates generated with k^j set to unity. For a fc^ variation of 1.6%, the 
variation of the fitted W mass in the EC is shown in Table 0. 

TABLE IV. Value of fc^ giving the minimum for different pdf's. 



MRS(A') CTEQ3M g CTEQ2M [g 


MRSD-' m 


0.975 0.98 0.985 


0.99 


TABLE V. Variation in fitted EC W mass due to a 1.6% 


variation in /c^. 


rriT fit Pt{&) ht 


PT{y) fit 


5Mw (MeV) 34 48 


25 



The comparison of the electron rj distribution between the data and the Monte Carlo 
provides a consistency check of the predicted W rapidity distribution, and hence of the pdf's. 



IV shows the values of fc„ at which 
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The measured k^j being consistent with unityf] sets an upper bound on the pdf uncertainty. 
While this constraint can potentially be much more powerful with higher statistics obtained 
in future data-taking, it is presently weaker than the uncertainty in the modern pdf's. 
Therefore we do not use this constraint to set our final W mass uncertainty due to pdf's. 
However, since our data used for this constraint are independent of the world data used to 
derive the pdf's, we have additional evidence that the uncertainty on the W mass due to 
the pdf's is not being underestimated. 



IX. BACKGROUNDS 
A. W — ^ TV — > evvu 

The decay W tu ^ evVv is topologically indistinguishable from W ev. It is 
included in the fast Monte Carlo simulation (Sec. 0). This decay is suppressed by the 
branching fraction for r evV (17.83 ± 0.08)% |T^, and by the lepton cuts. It accounts 
for 1% of the events in the W sample. 



B. Hadronic Background 



QCD processes can fake the signature of a W ^ eu decay if a hadronic jet fakes the 
electron signature and the transverse momentum balance is mismeasured. 
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FIG. 41. The spectra of a sample of events passing electron identification cuts (•) and a 
sample of events failing the cuts (o). 



^We have used A;^ = 1 in the mass analysis. 



44 



We estimate this background from the j^t spectrum of data events with an electromag- 
netic cluster. Electromagnetic clusters in events with low j^t are almost all due to jets. Some 
of these clusters satisfy our electron selection criteria and fake an electron. From the shape 
of the j^T spectrum for these events we determine how likely it is for these events to have 
sufficient to enter our W sample. 

We determine this shape by selecting isolated electromagnetic clusters that have > 
200 and the 4- variable likelihood A4 > 30. Nearly all electrons fail this cut, so that the 
remaining sample consists almost entirely of hadrons. We use data collected using a trigger 
without the requirement to study the efficiency of this cut for jets. If we normalize the 
background spectrum after correcting for residual electrons to the electron sample, we obtain 
an estimate of the hadronic background in an electron candidate sample. Figure ^ shows the 

spectra of both samples, normalized for -^t < 10 GeV. We find the hadronic background 
fraction of the total W sample after all cuts to be /had = (3.64 ± 0.78)%. The error receives 
contributions from the uncertainty in the relative normalization of the two samples at low ^t, 
the statistics of the failed electron sample, and the uncertainty in the residual contamination 
of the failed electron sample by true electrons. We fit the distributions of the background 
events with > 30 GeV to estimate the shape of the background contributions to the 
PT(e), Pt{^), and spectra (Fig. ^). We use the statistical error of the fits to estimate 
the uncertainty in the background shapes. 
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FIG. 42. Shapes of m^, pri^), and pt{i^) spectra from hadron ( ) and Z boson ( ) 

backgrounds with the proper relative normalization. 



C. Z ^ ee 

To estimate the fraction oi Z ^ ee events that satisfy the W boson event selection, 
we use a Monte Carlo sample of approximately 100,000 Z ^ ee events generated with 
the HERWIG program and a detector simulation based on geant. The boson px spectrum 
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generated by HERWIG agrees reasonably well with the calculation in Ref. ||29| and with our Z 
boson "px measurement |^|. Z ^ ee decays typically enter the W sample when one electron 



satisfies the W cuts and the second electron is lost or mismeasured, causing the event to 
have large ^t- 

An electron is most frequently mismeasured when it goes into the regions between the 
CC and one of the ECs, which are covered only by the hadronic section of the calorimeter. 
These electrons therefore cannot be identified, and their energy is measured in the hadronic 
calorimeter. Large is more likely for these events than when both electrons hit the EM 
calorimeters. 

We make the W and Z selection cuts on the Monte Carlo events, and normalize the 
number of events passing the W cuts to the number of W data events, scaled by the ra- 
tio of selected Z data and Monte Carlo events. We estimate the fraction of Z events in 
the W sample to be jz = (0.26 ± 0.02) %. The uncertainties quoted include systematic 
uncertainties in the matching of momentum scales between Monte Carlo and collider data. 
Figure ^ shows the distributions of Pt{g), Pt{i^)i and ttit for the Z events with one lost or 
mismeasured electron that satisfy the W selection. 



X. MASS FITS 

A. Maximum Likelihood Fitting Procedure 

We use a binned maximum likelihood fit to extract the W mass. Using the fast Monte 
Carlo program, we compute the m^, pr(e), and pri.^) spectra for 200 hypothesized values 
of the W mass between 79.7 and 81.7 GeV. For the spectra we use 250 MeV bins. The 
statistical precision of the spectra for the W mass fit corresponds to about 8 million W 
decays. When fitting the collider data spectra, we add the background contributions with 
the shapes and normalizations described in Sec. to the signal spectra. We normalize 
the spectra within the fit interval and interpret them as probability density functions to 
compute the likelihood 

N 

L{m) = X{pT{m), (42) 

i=l 

where Pi{m) is the probability density for bin i, assuming Mw = tu, and is the number 
of data entries in bin i. The product runs over all bins inside the fit interval. We fit 
— ln[L(m)] with a quadratic function of m. The value of m at which the function assumes 
its minimum is the fitted value of the W mass and the 68% confidence level interval is the 
interval in m for which — ln[L(m)] is within half a unit of its minimum. 



B. Electron pt Spectrum 



We fit the Pt{g) spectrum in the region 32 < pxie) < 50 GeV. The interval is chosen to 
span the Jacobian peak. The data points in Fig. ^ represent the pri^) spectrum from the 



46 



W sample. The solid line shows the sum of the simulated W signal and the estimated back- 
ground for the best fit, and the shaded region indicates the sum of the estimated hadronic 
and Z — > ee backgrounds. The maximum likelihood fit gives 

Mw = 80.547 ± 0.128 GeV (43) 



for the W mass. Figure H4| shows — ln(L(m) / Lq) for this fit, where Lq is an arbitrary number. 



. 700 




45 50 55 

p^Ce) (GeV) 

FIG. 43. Spectrum of pT(e) from the W data. The superimposed curve shows the maximum 
Ukelihood fit and the shaded region the estimated background. 

As a goodness-of-fit test, we divide the fit interval into 0.5 GeV bins, normalize the 
integral of the probability density function to the number of events in the fit interval, and 
compute = Y^iLiiVi ~ Pi?' IVi- The sum runs over ah bins, yi is the observed number 
of events in bin and is the integral of the normalized probability density function over 
bin %. The parent distribution is the distribution for N — 2 degrees of freedom. For the 
spectrum in Fig. ESIwe compute = 46. For 36 bins there is a 8% probability for > 46. 
Figure ^ shows the contributions Xi = iVi ~ Pi) I Jlh to x^ for the 36 bins in the fit interval. 



Figure ^ shows the sensitivity of the fitted mass value to the choice of fit interval. The 
points in the two plots indicate the observed deviation of the fitted mass from the value 
given in Eq. |4^. We expect some variation due to statistical fluctuations in the spectrum 
and systematic uncertainties in the probability density functions. We estimate the effect 
due to statistical fluctuations using Monte Carlo ensembles. We expect the fitted values to 
be inside the shaded regions indicated in the two plots with 68% probability. The dashed 
lines indicate the statistical error for the nominal fit. Figure ^ shows that the probability 
density function provides a good description of the observed spectrum. 
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FIG. 44. The likelihood function for the pT{e) ht. 
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FIG. 45. The x distribution for the fit to the pt{^) spectrum. 
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FIG. 46. Variation of the fitted mass with the prie) fit window limits. See text for details. 



C. Transverse Mass Spectrum 



The rriT spectrum is shown in Fig. 14^. The points are the observed spectrum, the solid 
line shows signal plus background for the best fit, and the shaded region indicates the 
estimated background contamination. We fit in the interval 65 < rriT < 90 GeV. Figure |8 
shows — ln(L(m)/Lo) for this fit where Lq is an arbitrary number. The best fit occurs for 

Mw = 80.757 ± 0.107 GeV. (44) 
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FIG. 47. Spectrum of mx from the W data. The superimposed curve shows the maximum 
likelihood fit and the shaded region shows the estimated background. 
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FIG. 48. The likelihood function for the niT fit. 



Figure ^ shows the deviations of the data from the fit. Summing over all bins in the 
fitting window, we get = 17 fo^' 25 bins. For 25 bins there is a 81% probability to obtain 



a larger value. Figure 50 shows the sensitivity of the fitted mass to the choice of fit interval. 
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FIG. 49. The x distribution for the fit to the rriT spectrum. 



D. Neutrino pt Spectrum 



Figure ^ shows the neutrino pt spectrum. The points are the observed spectrum, the 
solid line shows signal plus background for the best fit, and the shaded region indicates the 
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FIG. 50. Variation of the fitted mass with the m-r fit window limits. See text for details. 



estimated background contamination. We fit in the interval 32 < Pt(^) < 50 GeV. Figure p2 
shows — ln(L(m)/Lo) for this fit where Lq is an arbitrary number. The best fit occurs for 



w 



80.740 ±0.159 GeV. 



(45) 




Pt(v) (GeV) 

FIG. 51. Spectrum of pt(z^) from the W data. The superimposed curve shows the maximum 
likelihood fit and the shaded region shows the estimated background. 



Figure ^ shows the deviations of the data from the fit. Summing over all bins in the 
fitting window, we get = 37 for 36 bins. For 36 bins there is a 33% probability to obtain 



a larger value. Figure 54 shows the sensitivity of the fitted mass to the choice of fit interval. 



51 



o 

^ - 




80 80.5 81 81.5 



(GeV) 

FIG. 52. The likelihood function for the pt{v) fit. 
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FIG. 53. The x distribution for the fit to the pt{i^) spectrum. 
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FIG. 54. Variation of the fitted mass with the pt{i^) fit window limits. See text for details. 

XI. CONSISTENCY CHECKS 

A. North vs South Calorimeters 

Since the detector is north-south symmetric, we expect the measurements made with the 
north and south calorimeters separately to be consistent. We find 



ECN 
W 



ECN 
W 

fECN 



ECS 
W 
ECS 
W 



ECS 



88 ± 215 MeV (mr fit) 
-116 ± 258 MeV (p^ fit) 
107 ± 318 MeV {p^ fit) 



(46) 



where the uncertainty is statistical only. 



B. Time Dependence 

We divide the W boson data sample into five sequential calender time intervals such 
that the subsamples have equal number of events. We generate resolution functions for the 
luminosity distribution of these five subsamples. We fit the transverse mass and lepton 
spectra from the W samples in each time bin. The fitted masses are plotted in Fig. ^ 
where the time bins are labelled by run blocks. The errors shown are statistical only. We 
compute the with respect to the W mass fit to the entire data sample. The per degree 
of freedom (dof) for the Pt{g) fit is 7.0/4 and for the J9t(^) fit is 1.5/4. The rriT fit has a 
xVdof of 2.1/4. 

Since the luminosity was increasing with time throughout the run, the time slices corre- 
spond roughly to luminosity bins. 
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FIG. 55. The fitted W boson masses in bins of run blocks from the rriT, Pt{^)-, and pt{i^) fits. 
The solid line is the central value for the respective fit over the entire sample. The W fit statistical 
error for each subsample is shown. The average instantaneous luminosity in the bins is 4.2, 6.1, 
7.1, 9.3 and 10.1 respectively, in units of lO'^'^/cm^/s. 



C. Dependence on ut Cut 



We change the cuts on the recoil momentum ut and study how well the fast Monte Carlo 
simulation reproduces the variations in the spectra. We split the W sample into subsamples 
with M|j > GeV and u\\ < GeV, and fit the subsamples with corresponding Monte Carlo 
spectra generated with the same cuts. The difference in the fitted masses from the two 
subsamples corresponds to 0.3a, 0.8a and l.So" for the my, Pt{g), and Pt{^) fits respectively, 
based on the statistical uncertainty alone. Although there is significant variation among the 
shapes of the spectra for the different cuts, the fast Monte Carlo models them well. 



D. Dependence on Fiducial Cuts 

We fit the mrp spectrum from the W sample and the m{ee) spectrum from the Z sample 
for different pseudorapidity cuts on the electron direction. Keeping the upper |?7det(e)| cut 
fixed at 2.5, we vary the lower |77det(e)| cut from 1.5 to 1.7. Similarly, we vary the upper 
|^det(e)| cut from 2.0 to 2.5, keeping the lower |?]det(e)| cut fixed at 1.5. Figures |56|-p^ show 
the change in the W mass versus the ?7det(e) cut using the electron energy scale calibration 
from the corresponding Z sample. The shaded region indicates the statistical error. Within 
the uncertainties, the mass is independent of the ?7det(e) cut. 
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FIG. 56. The variation in the W mass from the pt{g) fit versus the 7ydet(e) cut. The shaded 
region is the expected statistical variation. 
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FIG. 57. The variation in the W mass from the mT fit versus the ?/dct(e) cut. The shaded 
region is the expected statistical variation. 
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FIG. 58. The variation in the W mass from the Pt{i^) fit versus the ?7det(e) cut. The shaded 
region is the expected statistical variation. 

E. Z Boson Transverse Mass Fits 

As a consistency check, we fit the transverse mass distribution of the Z — >■ ee events, re- 
constructed using each electron and the recoil. The measured energy of the second electron is 
ignored, both in the data and in the Monte Carlo used to obtain the templates. Each Z event 
is treated (twice) as a event, where the neutrino transverse momentum is recomputed us- 
ing the first electron and the recoil. One of the two electrons is required to be in the EC. The 
fitting range is 70 < < 90 GeV for the CC/EC events and 70 < my < 100 GeV for the 
EC/EC events. Figure |U| shows the results. The CC/EC fit yields Mz = 92.004±0.895(stat) 
GeV with xVdof = 7/9. The EC/EC fit yields Mz = 91.074 ± 0.299(stat) GeV with xV^of 
= 16/14. The average fitted mass is Mz = 91.167 ± 0.284(stat) GeV. The fits are good and 
the fitted masses are consistent with the input Z mass. 
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XII. UNCERTAINTIES IN THE MEASUREMENT 

Apart from the statistical error in the fitted W mass, uncertainties in the various inputs 
needed for the measurement lead to uncertainties in the final result. Some of these inputs 
are discrete (such as the choice of the parton distribution function set) and others are 
parameterized by continuous variables. For a different choice of pdf set, or a shift in the 
value of an input parameter by one standard deviation, the expected shift in the fitted 
W mass is computed by using the fast Monte Carlo to generate spectra with the changed 
parameter and fitting the spectra with the default templates. The expected shifts due to 
various input parameter uncertainties (given in Table [V|) or choice of pdf set are discussed 
in detail below, and are summarized in Tables |V11| and |V111| . The shifts in the fitted mass 



obtained from the different kinematic spectra may be in opposite directions, in which case 
they are indicated with opposite signs. 
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FIG. 59. Spectra of the Z boson transverse mass, from the CC/EC data (top) and the EC/EC 
data (bottom). The second electron in the Z boson decay is treated hke the neutrino in W boson 
decay. The superimposed curves show the maximum hkehhood fits and the shaded regions show 
the estimated backgrounds. The x^/dof between the data and the Monte Carlo are also shown. 
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TABLE VI. Errors on the parameters in the W mass analysis. The correlation coefficient 
between a^ec and Prec is —0.98; that between Srec and Omb is —0.60. 



parameter 


error 


parton luminosity (3 


0.001 GeV -1 


photon coalescing radius Rq 


7 cm 


W width 


59 MeV 


ECEM offset 6ec 


0.7 GeV 


ECEM scale qec 


0.00187 


FDC radial scale /3fdc 


U.UUU04 


FDC-EC radial scale /3ec 


0.0003 


ECEM constant term cec 


+U.UU6 
-0.01 


recoil response {oj-ec, Prec) 


(0.06, 0.02) 


recoil resolution (sreo Omb) 


(0.14 GeV^/^, 0.028) 




e (0.0,0.01) 


uy correction Au^^/5r]6(j) 


28 MeV 


U|| efficiency slope s 


0.0012 GeV'i 



Since the most important parameter, the EM energy scale, is measured by calibrating 
to the Z mass, we are measuring the ratio of the W and Z boson masses. There can be 
significant cancellation in uncertainties between the W and Z masses if their variation due 
to an input parameter change is very similar. For those parameters that affect the fitted Z 
mass. Tables [VI I| and [VIII| also show the expected shift in the fitted Z mass. The signed W 
and Z mass shifts are used to construct a covariance matrix between the various fitted W 
mass results, which is used to obtain the final W mass value and uncertainty; thus simple 
combination of the uncertainties in Tables |VII| and |VIII| is inappropriate. This is discussed 



in detail in Section XIII 



A. Statistical Uncertainties 



Tables |V11| and |V111| list the uncertainties in the W mass measurement due to the finite 



sizes of the W and Z samples used in the fits to the m^, Pxi^), Pri^^), and m{ee) spectra. The 
statistical uncertainty due to the finite Z sample propagates into the W mass measurement 
through the electron energy scale oec- 

Since the m-r, Pri^) and pr(i^) fits are performed using the same W data set, the results 
from the three fits are statistically correlated. The correlation coefficients between the 
respective statistical errors are calculated using Monte Carlo ensembles, and are shown in 
Table ^ 



B. W Boson Production and Decay Model 



58 



TABLE VII. Variation in the fitted Mw and Mz (in MeV) for the forward electron sample due 
to variation in the model input parameters by the respective uncertainties. 



Source 


5Mz 


6Mz 


6Mw 


6Mw 


5Mw 




(CC/EC) 


(EC/EC) 


{itlt) 


{p't) 


(p't) 


statistics 


124 


221 


107 


128 


159 


Pt{W) spectrum 






22 


37 


44 


MRSR2 |7| 






-11 


-21 


-43 


MRS(A') II 






-7 


-43 


-19 


CTEQ5M [||] 






14 


9 


-17 


CTEQ4M [||] 






1 


-21 


22 


CTEQ3M Q 






13 


30 


28 


parton 












luminosity /3 


8 


7 


9 


11 


18 


Ro 


10 


13 


9 


17 


12 


27 


5 


10 


5 


10 





W width 






10 


10 


10 


ECEM offset 


284 


421 


437 


433 


386 


ECEM scale 












variation 0.0025 


114 


228 


201 


201 


201 


CCEM scale 












variation 0.0008 


37 














FDC radial scale 


8 


36 


43 


37 


28 


FDC-EC radial scale 


10 


52 


57 


54 


48 


ECEM constant 












term cec 





U 


45 


on 
zy 


7o 


hadronic 












response 






11 


20 


-50 


hadronic 












resolution 






40 


4 


203 


correction 


20 


30 


18 


34 


-6 


efficiency 






4 


-22 


40 


background 












normalization 





11 


12 


15 


25 


background 












shape 





5 


16 


23 


78 
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TABLE VIII. Variation in the fitted Mw and Mz (in MeV) for the central electron sample 
due to variation in the model input parameters by the respective uncertainties. 



Source 


6Mz 


6Mz 






5Mv/ 




(CC/CC) 


(CC/EC) 


(mr) 


(Pt) 


(Pt) 


statistics 


75 


124 


70 


85 


105 


PtO^) spectrum 






10 


50 


25 


MRSR2 g7| 






5 


26 


3 


MRS (A') [Hi 






-5 


16 


-31 


CTEQ5M [|8| 






-8 


6 


-22 


CTE04M [H9II 






10 


11 


— 18 


CTEQ3M g 









64 


-9 


parton 












luminosity /3 


4 


8 


9 


11 


9 


Ro 


19 


10 


3 


6 





27 


10 


5 


3 


6 





W width 






10 


10 


10 


CC EM offset 


387 


467 


367 


359 


374 


CDC scale 


29 


33 


38 


40 


52 


uniformity 






10 


10 


10 


CCEM constant 












term cqc 






23 


14 


27 


hadronic 












response 






20 


16 


-46 


hadronic 












resolution 






25 


10 


90 


n|| correction 






15 


15 


20 


efficiency 






2 


-9 


20 


backgrounds 






10 


20 


20 
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TABLE IX. The statistical correlation coefficients obtained from Monte Carlo ensemble tests 
fitting the W boson mass for 260 samples of 11,089 events each. 





rriT 


correlation matrix 
VT{e) 




rriT 


1 


0.634 


0.601 


PT{e) 


0.634 


1 


0.149 




0.601 


0.149 
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1. Sources of Uncertainty 



Uncertainties in the W boson production and decay model arise from the following 
sources: the phenomenological parameters in the calculation of the priW) spectrum, the 
choice of parton distribution functions, radiative decays, and the W boson width. In the 
following we describe how we assess the size of the systematic uncertainties introduced by 



each of these. We summarize the size of the uncertainties in Tables VII and VIII 



2. W Boson pt Spectrum 
In Sec. |VII1| of Ref . M , we described our constraint on the W boson px spectrum. This 



constraint was obtained by studying the Z boson pt spectrum, which can be measured well 
using the two electrons in Z —>■ ee decays. For any chosen parton distribution function, the 
parameters of the theoretical model were tuned so that the predicted Z boson pt spectrum 
after simulating all detector effects agreed with the data. The precision with which the 
parameters could be tuned was limited by the statistical uncertainty and the uncertainty in 
the background. These parameter values were used to predict the W boson px spectrum. 
The uncertainties in the fitted W boson mass for the CC W sample due to the uncertainty 



in the W boson p^ spectrum were listed in Ref. [Q, and are reproduced in Table |VIII| . The 
corresponding uncertainty in the EC analysis is given in Table |V11[ The CC and EC W 
mass uncertainties from this source are assumed to be fully correlated. 



3. Parton Distribution Functions 

To quantify the W mass uncertainty due to variations in the input parton distribu- 
tion functions, we select the MRS(A'), MRSR2, CTEQ5M, CTEQ4M and CTEQ3M sets 
to compare to MRST. We select these sets because their predictions for the lepton charge 
asymmetry in W decays and the neutron-to-proton Drell-Yan ratio span the range of consis- 
tency with the measurements from CDF and E866 [^. These measurements constrain 



the ratio of u and d quark distributions which have the most influence on the W rapidity 
spectrum. 
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Using these parton distribution function sets as input to the fast Monte Carlo model, 
we generate rriT and lepton spectra. For each chosen parton distribution function set we 
use the appropriate W boson pt spectrum as used in our CC W mass analysis. We then fit 
the generated spectra in the same way as the spectra from collider data, i.e. using MRST 
parton distribution functions. Table |V11| lists the variation of the fitted EC W mass values 
relative to MRST. The CC and EC W mass uncertainty from this source is taken to be fully 
correlated, taking the relative signs of the mass shifts into account. 

We find that the combination of the CC and EC W boson mass measurements is less 
sensitive to pdf variations, than for the CC measurement alone. The pdf uncertainty on the 
CC measurement is 11 MeV. The pdf uncertainty on the CC+EC combined measurement is 
7 MeV. As expected, the larger combined rapidity coverage makes the observed transverse 
mass and transverse momentum distributions less sensitive to the longitudinal boost of the 
W boson. 



4- Parton Luminosity 

The uncertainty of 10^^ GeV^^ in the parton luminosity slope /3 (Sec. 0) translates into 
an uncertainty in the fitted W and Z boson masses. We estimate the sensitivity in the fitted 
W and Z masses by fitting Monte Carlo spectra generated with different values of f3. The 
uncertainty in f] is taken to be fully correlated between the CC and EC W mass analyses. 



5. Radiative Decays 

We assign an error to the modeling of radiative decays based on varying the detector 
parameter Rq (Sec. 0). Rq defines the maximum separation between the photon and electron 
directions above which the photon energy is not included in the electron shower. In general, 
radiation shifts the fitted mass down for the transverse mass and electron fits, because for a 
fraction of the events the photon energy is subtracted from the electron. Hence increasing Rq 
decreases the radiative shift. Both the fitted W and Z masses depend on Rq. To estimate the 
systematic error, we fit Monte Carlo spectra generated with different values of Rq. geant 
detector simulations show that, for an Rq variation of ±7 cm, the electron-photon cluster 
overlap changes to give the maximum variation in the electron identification efficiency. The 
changes in the mass fits when varying Rq by ±7 cm are listed in Table |V11| . 

There are also theoretical uncertainties in the radiative decay calculation. Initial state 
QED radiation is not included in the calculation of Ref. PDI. However, initial state radiation 



does not affect the kinematic distributions used to fit the mass in the final state. We 
studied the effect of QED radiation off the initial state quarks on the parton luminosity by 
computing the parton luminosity including and excluding QED radiative effects on the quark 
momentum distribution. The change in the parton luminosity slope parameter was less then 
half of the quoted uncertainty on the parameter, which was dominated by acceptance effects. 
The calculation of Ref. |^ includes only processes in which a single photon is radiated. 



We use the code provided by the authors of Ref. |Q to estimate the shift introduced in the 
measured W and Z masses by neglecting two-photon emission. The estimated shifts in the 
W and Z fitted masses due to two-photon radiation are shown in Table |V11|. Since this effect 
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is an order of magnitude smaller than the statistical uncertainty in our measurement we do 
not correct for it, but add it in quadrature to the uncertainty due to radiative corrections. 
The uncertainty in the radiative correction is taken to be fully correlated between the CC 
and EC W mass analyses. 



6. W Boson Width 

The uncertainty on the fitted W mass corresponds to the uncertainty in the measured 
value of the W boson width Tw = 2.062±0.059 GeV H- We take this uncertainty to be 
fully correlated between the CC and EC W mass analyses. 



Our recent measurement of the W width [54| considerably improves the precision of F^y 



and would reduce the W mass uncertainty from this source. However, since this is already 
a small source of uncertainty, the impact on the total W mass uncertainty is small. 



C. Detector Model Parameters 



The uncertainties on the parameters of the detector model determined in Sees. |VI|- |VI1 



translate into uncertainties in the W mass measurement. We study the sensitivity of the W 
mass measurement to the values of the parameters by fitting the data with spectra generated 
by the fast Monte Carlo with input parameters modified by ±1 standard deviation. 

Table IVllI lists the variation in the measured EC W mass due to variation in the indi- 



vidual parameters. For each item the uncertainty is determined with a typical Monte Carlo 
statistical error of 5 MeV. To achieve this precision, 10-20 million W —>■ eu decays are 
simulated for each item. 

The residual calorimeter nonlinearity is parametrized by the offset 5ec- The electron 
momentum resolution is parametrized by cem- The electron angle calibration includes the 
effects of the parameters Pfdc and Pec, discussed in Section |V|. The recoil response is 
parameterized by ctrec and P^ec- The recoil resolution is parameterized by s^ec and Omb- 
Electron removal refers to the bias Au\\ introduced in the u\\ measurement by the removal 
of the cells occupied by the electron. Selection bias refers to the m efficiency. 



D. Backgrounds 

We determine the sensitivity of the fit results to the assumed background normalizations 
and shapes by repeating the fits to the data with background shapes and normalizations 
modified by ±1 standard deviation. Table [VI I| lists the uncertainties introduced in the EC 
W boson mass measurement. 



XIII. COMBINED EC AND CC W BOSON MASS ERROR ANALYSIS 

The measurement of the W mass requires the knowledge of many parameters in our model 
of the W production, decay and detector response. These parameters are constrained by 
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measurements, and in some cases by theoretical input. The W mass error analysis involves 
the propagation of the measurement or theoretical uncertainties to the error matrix on the 
parameters, which is then propagated further to the error matrix on the CC and EC W 
mass measurements. The error matrix allows us to combine the fitted W mass values using 
the different data samples and techniques into a single value with a combined error. 

Wc identify the following parameters of relevance to the W mass measurements in the 
EC and CC: 

W mass statistical errors Sou cc ^-nd Sujec 

EM scales acc and a^c 

EM offset parameters See and Sec 

FDC scale /3fdc and FDC-EC relative scale /3ec 

CDC scale Pcdc 

EM resolutions (constant terms) ccc and cec 

recoil response Orec representing jointly the response parameters CKj-ec and /9rec 

recoil resolution g^cc representing jointly the hadronic sampling term Srec and the effects 
of the underlying event a^ah 

backgrounds bcc and 6ec 

It 1 1 corrections ucc and uec 

M|| efficiencies Eqc and Eec 

radiative corrections as a function of the photon coalescing radius Rq 
parton luminosity /3 
theoretical modeling t 

We take the EM scales, EM offsets, angular scales, u\\ corrections, parton luminosity and 
the radiative correction to be a set of parameters that jointly determine the measured W 
and Z masses. We also take the EM resolution parameters as a correlated set. We take 
the CC and EC backgrounds and u\\ efficiencies to be uncorrelated. The recoil modeUing 
and the theoretical modelling (including pdf 's, Pt{W) spectrum, parton luminosity, radiative 
corrections and W width) are treated as being common between the CC and the EC analyses. 
For all correlated parameters the sign of the W mass correlation is determined by the relative 
sign of the mass shifts. 

The following measurements provide information on the values of these parameters 

• The Z mass measurements M^*"^*"*", M^*"^^*" and M^*"^^*" 



• FDC radial calibration ^pDC and FDC-EC relative radial cahbration ^ec 
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• CDC z calibration ^cdc 

• CC and EC EM offset measurements occ and oec 

• Gaussian widtli fitted to Z boson peak cr^'"^'"^, o^^^^ and o^^^'^ 

• pt balance in Z events 

• width of balance in Z events 

• measurements of My correction and My efficiency 

• constraints on theoretical model (boson from D0 data, W width from world data 
including D0 data, and pdf 's and parton luminosity from world data) 

We express the variations on the various calibration quantities (such as Z mass, EM 
offset, and angular scales, collectively referred to as C) and the Z width measurements as a 
linear combination of the variations on the parameters 

5az = A, 5cEM (47) 

where 

5C = (5mJ^/^^, 5mJ^/^^, 5Mz^^^^, 5^fdc, ^^ec, 

Sp= (5«CC, ^ttEC, 5/3fDC, 5/3eC, 5/3cDC, 

S5cc,SSec,SRo,Succ,Suec,S(3) (48) 

r — * / r CC/CC r CC/EC r EC/EC\ 

Saz = [daz ,oaz ' ,daz ), 
(^cem = ((^ccc, (^cec)- (49) 

The A matrices contain the partial derivatives of the observables with respect to the pa- 
rameters. 

Similarly, the variations on the W mass are related linearly to the parameter variations 

5Mw = Aiy 5p 

~l~ Aj-Qcoil scale '^'^rec 

+ A 

recoil resolution 

— * 

^background 

+ Au 5u 
+ A, 8e 

+ Atheory 

+ 5u (50) 



and 
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where 5Mw = {5M^^,5M^^). 

Knowing the components of 6C and Saz, we compute the covariance matrix for the 
parameters in p and cem- Since there are more measurements than parameters, we use 
the generahzed least squares fitting procedure for this purpose. We then propagate the 
parameter covariance matrices into the covariance matrix for the CC and EC W mass 
measurements using equation by identifying the covariance matrix with the expected 
value of 6Mw{SMw)'^ , where T indicates the transpose. The various contributions to 6Mw 
are independent, hence they contribute additively to the total covariance matrix. 

The CC W mass measurements were obtained using the MRS (A') parton distribution 
functions. We adjust these measurements by the estimated shifts (see Table |V111| ) when 
using the MRST parton distribution functions. Thus we use the following W mass values 
extracted from the CC data to combine with our EC measurements: 

M^^ = 80.443 GeV (m^ fit) 
M^^ = 80.459 GeV (prie) fit) 

M^^ = 80.401 GeV (pt(i^) fit) (51) 

The combined W mass Mw for a set of n W mass measurements rrii and their covariance 
matrix V is given by 

n n 

Mw = iJ2 H^,m^)/{Y. H,-), (52) 

i,j=l i,j=l 

where H = and i,j run over the W mass measurements being combined. The combined 
error is given by 

a{M^) = {f:H,,)-'/\ (53) 
and the for the combination is given by 

n 

X^= J2 K - Mw) Hi, (m, - Mw). (54) 



XIV. RESULTS 

We use the covariance matrix described above to obtain the total uncertainty on the 
EC W mass measurements and to combine our CC and EC measurements. We obtain the 
following results for the transverse mass fit 

M^^ = 80.757 ± 0.107(stat) ± 0.204(syst) GeV 

= 80.757 ± 0.230 GeV (55) 



and 
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Mw = 80.504 ± 0.097 GeV (CC and EC combined). (56) 

The foi' the CC+EC ttit combination is 1.5 for one degree of freedom, with a probabihty 
of 23%. 

Similarly, for the pT(e) fit we obtain 

M^^ = 80.547 ± 0.128(stat) ± 0.203(syst) GeV 

= 80.547 ± 0.240 GeV (57) 

and 

Mw = 80.480 ± 0.126 GeV (CC and EC combined). (58) 

The fo^^ the CC+EC pr(e) combination is 0.1 with a probability of 74%. 
For the Pt{i^) fit we obtain 

M^^ = 80.740 ± 0.159(stat) ± 0.310(syst) GeV 

= 80.740 ± 0.348 GeV (59) 

and 

Mw = 80.436 ± 0.171 GeV (CC and EC combined). (60) 

The for the CC+EC Pt{i^) combination is 1.0 with a probability of 32%. 

The combination of the ttit, Pr(e) and Pt{^) fit values for the EC give the combined EC 
W mass result 

Mw = 80.691 ± 0.227 GeV. (61) 

The xVdof is 4.0/2, with a probability of 14%. 

We combine all six measurements (CC and EC fits with the three techniques) to obtain 
the combined 1994-1995 measurement 

Mw = 80.498 ± 0.095 GeV. (62) 

The x^/dof is 5.1/5, with a probability of 41%. The consistency of the six results indicates 
that we understand the ingredients of our model and their uncertainties. Including the 
measurement from the 1992-1993 data gives the 1992-1995 data measurement: 

Mw = 80.482 ± 0.091 GeV. (63) 

Table ^ lists the D0 W mass measurement uncertainties from the 1994-1995 end 
calorimeter data alone and the combined 1994-1995 central and end calorimeter data. 

The D0 measurement is in good agreement with other measurements and is more precise 
than previously published results. Table ^ lists previously published measurements with 
uncertainties below 500 MeV, except previous D0 measurements which are subsumed into 
this measurement. A global fit to all electroweak measurements excluding the direct W 
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TABLE X. W mass uncertainties (in MeV) in the EC measurement and the combined CC+EC 
measurement from the 1994-1995 data. 
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FIG. 60. A comparison of this measurement with previously published W boson mass mea- 
surements (Table pC| ). The shaded region indicates the predicted W boson mass value from global 
fits to all electroweak data except the W mass measurements |11|. 
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TABLE XI. Previously published measurements of the W boson mass. 



measurement Mw (GeV) reference 

CDF 90 79.910±0.390 H 

UA2 92 80.360±0.370 

CDF 95 80.410±0.180 

L3 99 80.610±0.150 

ALEPH 99 80.423±0.124 0] 

OPAL 99 80.380±0.130 

DELPHI 99 80.270±0.145 

D0 99 combined (this result) 80.482±0.091 



mass measurements predicts Mw = 80.367 ± 0.029 GeV ||TTI]. Figure |63 gives a graphical 
representation of these data. 

We evaluate the radiative corrections Atew, defined in Eq. Our measurement of Mw 
from Eq. |63| leads to 



Avew = -0.0322 ± 0.0059, 



(64) 



5.5 standard deviations from the tree level value, demonstrating the need for higher-order 
electroweak loop corrections. In Fig. Bl] we compare the measured W boson and top quark 



masses [20| from D0 with the values predicted by the standard model for a range of Higgs 
mass values [Q. Also shown is the prediction from the calculation in Ref. for a model 
involving supersymmetric particles assuming the chargino, Higgs, and left-handed selectron 
masses are greater than 90 GeV. The measured values are in agreement with the prediction 
of the standard model, and in even better agreement with a supersymmetric extension of 
the standard model. 
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FIG. 61. A comparison of the W boson and top quark mass measurements by the D0 collab- 



oration with the standard model predictions for different Higgs boson masses |56|. The width of 
the bands for each Higgs boson mass value indicates the uncertainty due to the error in a(M^). 
Also shown is the range allowed by the MSSM [p2l. 
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